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ABSTRAGT 

Using the Hubble Space Telescope (i75'T)/WFC3 and Advanced Camera for Surveys multi-band 
imaging data taken in CANDELS and 3D-HST, we study the general properties and diversity of the 
progenitors of the Milky Way (MWs) and local massive galaxy (MGs) at 0.5 < z < 3.0, based on a 
constant cumulative number density analysis. After careful data reduction and stacking analysis, we 
conduct a radially resolved pixel spectral energy distribution fitting to obtain the radial distributions 
of the stellar mass and rest-frame colors. The stellar mass of MWs increases in self-similar way, 
irrespective of the radial distance, while that of MGs grows in an inside-out way where they obtain 
^ 75% of the total mass at outer (> 2.5 kpc) radius since z ^ 2. Although the radial mass profiles 
evolve in distinct ways, the formation and quenching of the central dense region (or bulge) ahead of 
the outer disk formation are found to be common for both systems. The sudden reddening of the 
bulge at 2 : ^ 1.6 and ^ ^ 2.4 for MWs and MGs, respectively, suggests the formation of the bulge 
and would give a clue to the different gas accretion histories and quenching. A new approach to 
evaluate the morphological diversity is conducted by using the average surface density profile and its 
dispersion. The variety of the radial mass profiles for MGs peaks at higher redshift {z > 2.8) and then 
rapidly converges to a more uniform shape at 2 ; < 1.5, while that for MWs remains in the outer region 
over the redshift. Compared with the observed star-formation rates and color profiles, the evolution 
of variety is consistently explained by the star-formation activities. 

Subject headings: galaxies: high-redshift — galaxies: evolution — galaxies: formation — galaxies: 
general — galaxies: structure — Galaxy: evolution 


1. INTRODUCTION 

The nature of high redshift (high-z) galaxies has 
been enthusiastically investigated for a decade. Thanks 
to the recent progress of observation studies with the 
Hubble Space Telescope (HST) and lOm-class ground- 
based telescopes, we have witnessed the clear and de¬ 
tailed structures of galaxies in the early Universe. Galax¬ 
ies in such an epoch are known to have quite differ¬ 
ent properties from local ones. One example is the so- 
called “red nugget^ ” which is passively evolving massive 

10^^ Mq) galaxy with small radius ^ 1 kpc (Daddi 
et al. 2007; Trujillo et al. 2007; van Dokkum et al. 2008; 
Damjanov 2009; Nelson et al. 2014). After a rapid for¬ 
mation phase at 2 : > 2, they are considered to increase 
stellar masses and sizes at relatively moderate speeds 
by accretions of satellite galaxies at the outer parts to 
evolve into massive elliptical galaxies in the local Uni¬ 
verse (Naab et al. 2009; Oser et al. 2010; Trujillo et 
al. 2012). 

On the other hand, the formation and evolution histo¬ 
ries of the bulge-disk structure of disk galaxies, including 
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the Milky Way (MW), are still matter of much debate; 
some studies explain them by in-situ star-formation with 
constant and moderate gas accretion from halo to cen¬ 
tral region (e.g., Kormendy & Kennicutt 2004; Dekel et 
al. 2009), while others use violent physical mechanisms, 
such as bar instability and giant clump migration (e.g., 
Noguchi 1998, 1999; Debattista et al. 2006). Although 
more elaborate studies, especially by means of galactic 
archaeology, have been conducted based on the chemical 
abundance and kinematics of resolved stars, the details 
of the formation and evolution history of the MW are 
still unknown (e.g., Lee et al. 2011; Bovy et al. 2012; 
Bergemann et al. 2014). 

Recently, high- 2 : progenitor studies have shed light on 
the early evolution of the local galaxies. By selecting 
progenitors that become, for example, massive elliptical 
galaxies in the local Universe, and observing their snap¬ 
shots over a certain redshift range, recent studies have 
investigated their evolution histories. Although there 
are a number of novel matching schemes to pick up pro¬ 
genitor galaxies (Ichikawa et al. 2007; Conroy & Wech- 
sler 2009; Guo et al. 2010; Tojeiro & Percival 2010), an 
example is based on a constant cumulative number den¬ 
sity, which ranks galaxies by physical quantities, such 
as stellar mass. The application of the analysis with 
constant cumulative number density to the MW progen¬ 
itors (MWs) at high- 2 : was conducted by van Dokkum 
et al. ( 2013 ) with HST imaging data. By stacking the 
i7-band images they derived average stellar mass pro¬ 
files and showed that MWs had evolved in a self-similar 
way, where the mass growth rates at the inner and outer 
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radii were comparable (see also Patel et al. 2013b). Com¬ 
parison with those of local massive galaxies (MGs) (van 
Dokkum et al. 2010; Patel et al. 2013a), which evolve in 
inside-out way, is also worth being mentioned (see also 
recent work by Papovich et al. 2014). However, the phys¬ 
ical interpretation for the self-similar evolution was not 
pursued in past studies, whereas the inside-out evolu¬ 
tion scenario was enthusiastically discussed (e.g., Oser et 
al. 2010; see also discussion in this study). 

In order to avoid a significant contribution from young 
stars, near-infrared light is favorably used to derive 
galaxy stellar mass profiles (Szomoru et al. 2013). Al¬ 
though H-hand images (or longer wavelength, Kauff- 
mann & Chariot 1998) are often applied to obtain 
the stellar mass profiles, more accurate ones could be 
achieved when the spatial variations of mass-to-light ra¬ 
tio, M/L, are reconstructed. Moreover, with light pro¬ 
files we could not know what physical mechanism(s) 
works in the galaxy evolution, whereas the studies of the 
local Universe seek for more specific evidences corrob¬ 
orated by, for example, mapping age distributions and 
chemical abundance of resolved stars. Multi-object inte¬ 
gral field units (IFUs) are also expected to achieve such 
scientific goals for galaxy evolution studies. However, 
such spectroscopic observations are still expensive, espe¬ 
cially for high- 2 : moderate mass galaxies of 10 ® M©. 
Probing the general properties of star-formation and ac¬ 
cumulation of stellar mass for such less massive galaxies 
still remains quite difficult. 

The only alternative is the spatially resolved, pixel-by- 
pixel spectral energy distribution (SED) modeling, which 
we utilize with some modifications. The pixel-by-pixel 
SED fit is not frequently used, because, for example, 
there are a limited number of available filter bands and 
a lack of spatial resolution with high signal-to-noise ra¬ 
tio (S/N). The first attempt of spatially resolved SED 
fitting for high-^ star-forming and clumpy galaxies was 
made by Wuyts et al. (2012) with the if5'T/Advanced 
Camera for Surveys (ACS) and WEC3 data (see also Zi- 
betti et al. 2009 for local galaxies). Including follow-ups 
(e.g., Guo et al. 2012), these studies showed that the 
star-forming region of high-^ galaxies was extended to 
outer radii, which is consistent with the recent obser¬ 
vations with lEU (Eorster Schreiber et al. 2014; Genzel 
et al. 2014). We note, however, that most of the tar¬ 
gets in the studies were massive (M^ > IO^^Mq) and 
actively star-forming galaxies at the moderate redshift 
range (0.7 < z < 1.5). Eurthermore, the galaxy images 
are binned non-uniformly to obtain high S/Ns. Such a 
binning method tends to blur morphological details. 

In the context, we make use of the constant cumu¬ 
lative number density method to select the progenitors 
of the MW and local massive galaxy (MWs and MGs, 
respectively). With highly resolved and deep multi¬ 
wavelength (optical to near-infrared) data provided by 
HST, we investigate the general properties and the di¬ 
versity of the progenitor galaxies, such as SEDs, and re¬ 
solved stellar mass and color profiles over the redshift of 
0.5 < 2 : < 3, or, from the cosmic high-noon to dusk 
(Madau et al. 1996; Heavens et al. 2004; Hopkins et 
al. 2006). Stacking the sample in each redshift bin, with 
correction for position angle and axis ratio, allows us to 
obtain the average light profiles of both populations with 


sufficient S/N for each pixel. With the multi-band av¬ 
erage profiles, we derive radially resolved SEDs, which 
give a clue about to the radial properties and evolution 
histories of stellar mass and rest-frame colors. Another 
novelty of this study is focusing on the diversity of galaxy 
morphology. By making full use of the derived average 
profile and individual profiles, we evaluate the morpho¬ 
logical varieties of the progenitors at each cosmological 
epoch for the first time. These varieties are then inves¬ 
tigated in conjunction with stellar mass, star-formation 
rate (SER) and rest-frame colors. 

The paper is organized as follows. In Section 2 we 
describe the data and catalog, as well as the selection 
method for the progenitors. The details of the resolved 
SED fitting and relevant reduction methods are followed 
in Section 3. The results and discussion are given in 
Sections 4 and 5. The paper concludes with a summary 
in Section 6. Throughout the paper, we assume = 
0.3, Da = 0.7 and Hq = 70 kms“^Mpc“^ for cosmologi¬ 
cal parameters and Chabrier (2003) initial mass function 
(IME) for SED models. We use the AB magnitude sys¬ 
tem (Oke & Gunn 1983; Eukugita et al. 1996). 

2. DATA 

2.1. Sample Selection from the 3D-HST Catalog 

The Gosmic Assembly Near-infrared Deep Extragalac- 
tic Legacy Survey (CANDELS; Grogin et al. 2011; 
Koekemoer et al. 2011) serves the multi-band optical 
to near-infrared imaging data of 2-10 orbits per point¬ 
ing with the HST for five well-known fields (GOODS- 
north and south, AEGIS, COSMOS and UDS). The to¬ 
tal survey area reaches ^900 arcmin^. These deep and 
highly resolved imaging data enable us to obtain suffi¬ 
cient S/Ns for each pixel, even for high- 2 ; (> 2) galax¬ 
ies, while eight-filter bands over optical to near-infrared 
wavelengths promise reliable stellar population synthe¬ 
ses. We make use of a publicly available galaxy catalog 
by 3D-HST team (v.4.1; Skelton et al. 2014). Although 
the full description of the catalog is omitted, a brief sum¬ 
mary is given here. The photometric redshift was ob¬ 
tained with the EAZY code (Brammer et al. 2008), which 
showed good consistency with previously obtained spec¬ 
troscopic redshifts, A^/(l-h 2 :) = 0.003 — 0.005, for galax¬ 
ies with H < 23 (Skelton et al. 2014). The stellar mass 
was obtained with the EAST code (Kriek et al. 2009), 
assuming a Chabrier (2003) IME. The best-fit SED was 
obtained with the HST imaging data, Spitzer/IRAC 
mid-infrared data, and ground-based imaging data at the 
wavelengths from ultraviolet U-band to near-infrared K- 
band. The SER was derived from the rest-frame UV 
flux of the best-fit SED template and mid-/far-infrared 
data by Spitzer/MlVS. In addition, medium-band imag¬ 
ing data are available for some regions, which increases 
the reliability of SED fitting. 

Erom the 3D-HST catalog, we select the sample based 
on the spectroscopic redshift if available and photomet¬ 
ric redshift for others, and stellar mass (M*). Eurther 
criteria for reliable photometry are set with star-flag=0 
and usejphot = 1, which limit the sample to the galaxies 
with H < 23 mag in E160W-band. Summing up image 
counts for higher redshift galaxies (therefore, with lower 
S/N) is an arduous and precarious task because wrong 
sky subtractions could lead to significantly biased results. 
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Fig. 1.— Progenitors of the Milky Way (MWs) and local mas¬ 
sive galaxy (MGs) extracted by Eqs. (1) and (2) (dashed and solid 
lines) shown with blue dots and red circles, respectively. The lim¬ 
its of 90 and 75% sample completeness for star-forming (SFGs, 
blue) and quiescent galaxies (QGs, red) are indicated on the or¬ 
dinate. A randomly selected sample for each population is de¬ 
picted in each redshift bin, composited with F814W (blue), F125W 
(green) and F160W (red) of the HST imaging data. The stel¬ 
lar mass and star-formation rate (SFR) of each sample, which 
are taken from the 3D-HST catalog, are shown in the inset in 
unit of logM© and logMQyr”^. All the samples are available at 
http: //www.astr.tohoku.ac.jp/~mtakahiro/moril5apj/mori_15_rg 
b.pdf. We follow MWs and MGs to 2 ; ~ 3, where they had (only) 
~ 7% and 29% of their present day masses (~ 4.6 x 10^® Mq and 
1.5 X 10^^Mq), respectively. 


especially for galaxies at the edge of observation fields, at 
nearby bright stars, in satellite trails, and in ^ 50 pixel 
circular dead pixels on the WFC3/IR detector (referred 
to as ’’Death Star” by the 3D-HST team). We visually 
exclude those erroneous galaxies. 

The selected progenitor galaxies are shown in Fig. 1, 
with a randomly selected galaxy at each sample bin. 
The sample is 90-95% complete jfor galaxies with > 
IO^^Mq out to 2 ; = 2.5 (van der Wei et al. 2014). It is 
noted that the adopted criteria are sufficient to fully in¬ 
clude the blue star-forming galaxies with ^ 10 ^ Mq 
at 2 : ^ 3, while some of red quiescent counterparts might 
be missed. However, since blue galaxies are dominant 
in numbers at higher redshift when the constant cumu¬ 
lative number density method is applied (van Dokkum 
et al. 2013), the sample completeness is large enough for 
our purpose. We show 90 and 75% sample complete¬ 
ness limits for star-forming (SFGs, blue) and quiescent 
galaxies (QGs, red), which are derived in van der Wei et 
al. (2014), in the figure. 

2 . 2 . Selection of the MW and MG Progenitors 

We follow the constant cumulative number density cri¬ 
terion derived by van Dokkum et al. (2013) and Patel et 
al. (2013a) to select, at each redshift, galaxies that are 
supposed to become the MW-like and massive galaxies 
at 2 ; ^ 0. By using the constant number density for each 
redshift, van Dokkum et al. (2013) derived the stellar 
mass growth equation as a function of redshift for MWs 
(n ^ 1.1 X 10“^ Mpc“^), 

log(M*/Mo) = 10.66 - 0.045Z - 0.13^^ (1) 


and Patel et al. (2013a) for MGs (n ^ 1.4X 10 Mpc ^), 

log(M*/M 0 ) = 11.19 - 0.068Z - 0.04z^ (2) 

where we convert the stellar mass with Kroupa (2001) 
IMF (M^^k) into one with Chabrier IMF (M^^c) through 
log(M*^c/^o) = log(M*^K/M^ 0 ) - 0.04 (Gimatti et 

al. 2008). The equations above are both derived with 
the mass functions of Marchesini et al. (2009) for high- 2 : 
galaxies, but at 2 : ^ 0 the former adopts the one of Mous- 
takas et al. (2013), while the latter of Gole et al. ( 2001 ). 
The difference between the two mass functions is negli¬ 
gible for the stellar mass range in the present study, and 
does not affect the sample selection. According to the 
equations, expected stellar masses are log(M>^/M 0 ) ^ 
9.51, 9.79, 10.04, 10.25, 10.40, 10.49, and 10.57 for MWs 
and log(M*/M 0 ) - 10.68, 10.79, 10.89, 10.96, 11.00, 
11.07, and 11.10 for MGs, respectively, at {z) ~ 2.8, 2.4, 
2.0, 1.6, 1 . 2 , 1.0, and 0.7. ^ Although the redshift bins 
have different numbers of sample (all bins include > 80 
for MWs and > 10 for MGs), the sample numbers are 
large enough to gain high S/N and for the following sta¬ 
tistical discussion. The numbers of galaxies used for the 
stacking are summarized in Table 1. It should be noted 
that there is difference in the sample numbers of F125W 
and F160W, both of which are thought to have similar 
observation depth and field coverage. By visually check¬ 
ing the imaging data, we find the difference originates 
from the satellite trails and bad pixels on the detector in 
either filter. In addition, some galaxies are hidden under 
the sky background at shorter wavelengths, due to the 
“LC-correction,” especially at higher redshift. We exclude 
such galaxies from the final sample. 

Since galaxies could change their ranks for several rea¬ 
sons (e.g., major merging and quenching), we cut the 
sample with a range of 0.2 dex in the stellar mass (see 
Leja et al. 2013). Although some of previous studies 
adopted different number densities for the progenitor se¬ 
lection (Leja et al. 2013; Ownsworth et al. 2014), we ver¬ 
ified that adopting such number densities do not change 
our results. The sample in this study is selected with 
box of ^ and M^, as well as previous studies, rather than 
in curved box where stellar mass criteria evolve continu¬ 
ously with redshift. Although the sample selection with 
box could bias the sample if the box size is large, in 
the present study we find no significant difference in the 
median values of the parameters (redshift, stellar mass, 
and SFR) between the two selection criteria. The me¬ 
dian values for M>^, 2 :, and SFR are also summarized in 
Table 2. 

3. METHODOLOGY 

Here we introduce the method of radially resolved SED 
fit, or radial SED, to derive the stellar mass and rest- 
frame colors for each pixel of the galaxy radial profiles. 
We also visually summarize the contents in this Section 
in Eig. 2. 

3.1. Data Assessment 

We use the imaging data of WEC3/IR (science, weight, 
and exposure images) reduced by the 3D-HST team 

® All the galaxies used for stacking are exhibited at 
http: / / www.astr.tohoku.ac.jp/~mtakahiro/moril5apj / mori_15_rg 
b.pdf. 
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(Skelton et al. 2014). For the data taken with ACS, how¬ 
ever, we have no access to the exposure maps made by 
the 3D-HST team. Therefore, for ACS imaging data we 
assume the exposure time for each, based on the sched¬ 
uled observation time (A. M. Koekemoer 2014, private 
communication). 

The error estimate of the flux is one of the most im¬ 
portant processes for deriving SEDs. Calculation for the 
ith pixel error, is made by the following equation, 

Ei = + (^ 0 /- (3) 

The first term of the right side, is the ith pixel 

value of the variance map (or inverse weight map, 
Avht.fits)^ in units of (e“/s)^, defined as 

_2 _ {Di ^ Bi) ^ Niaread^ ^ 

where Di is the total accumulated dark current signals 
and Bi is the total accumulated background level during 
the exposure on the ith pixel, both in unit of e“. Ni is the 
total number of exposures, cTread is the readout noise, and 
ti is the total exposure time. We note that the variance 
map does not include the Poisson noise of the objects, but 
only background (sky, dark current, and readout) noise, 
including the pixel correlation noise during the drizzle 
task (see Casertano et al. 2000 for details). Therefore, 
we need to add the Poisson noise for the object, (Jo^i^ in 
unit of (e“/s^), as defined by 

cTo/ = Fi/tl, (5) 

where Fi denotes the total photon count (e“) on the ith 
pixel of the reduced imaging data. 

The science images are then convolved with point- 
spread function (PSF) to be blurred to the FWHM of 
F160W (^ 0."18), while the error maps are convolved 
in quadrature. The convolution is conducted with con¬ 
volution kernels generated by the “CLEAN” algorithm 
(Hogbom 1974), where PSF images provided by the 3D- 
HST team (median stacked stars) are used. Examples 
of PSF profiles convolved to match the F160W PSF are 
shown in Fig. 3a, where we see that all the profiles are in 
good agreement, within ^ 2% difference, which is much 
smaller than the photometric error (see also appendix). 
At last we have two convolved image data, one for the 
science map (where the ith pixel has Fi) and the other 
for the error map (F^), from which we calculate the S/N 
of the ith pixel as Fi jEi, for each filter band. 

3.2. Conversion of Two-dimensional (2-D) Image to 
Onedimensional (1-D) Radial Light Profile 

Prior to stacking the galaxies in redshift bins, we con¬ 
vert the 2-D imaging data into 1-D radial profiles (here¬ 
after, ID-conversion) to obtain the average radial pro¬ 
file. We first extract the 2-D image of 500x500 pixel, 
or ^250x250 kpc at z ^ 2, from the original mosaic 
images to subtract the local sky background. Since the 
calculated sky background calculated would sometimes 
be overestimated by bright neighboring objects (e.g., 
HauBler et al. 2007), the following procedures are con¬ 
ducted. First, we aggressively mask all objects detected 
with SExtractor (Bertin & Arnouts 1996) and the pa¬ 
rameters in Table 3. The masked pixels are not included 


in the calculation of sky. We then estimate the local sky 
value from the median of the masked stamp. Although 
wrong sky subtractions would give systematically biased 
results, especially for the ground-based near-infrared ob¬ 
servations, the sky subtraction for the imaging data with 
HST would not affect the results in this study because 
sky is negligibly fainter than the stacked profiles. It is 
noted that applying the ID-conversion after stacking 2-D 
images could yield false spatial asymmetricities. 

In the 1-D conversion, we correct the effect of galaxy 
inclination. We obtain the projected axis ratio and apply 
the following equation (see also Bertin & Arnouts 1996) 
to convert 2-D imaging data into the inclination cor¬ 
rected 1-D profile, 

^^(^5 q) = (cos^ 0 + sin^ 0lq‘^){x — xo)‘^ 

+{siB^0 + cos^d/q^){y-yo)'^ ( 6 ) 

+2cos6»(l - l/q'^)ix - xo){y - yo), 

where 0 and q are position angle and axis ratio derived 
with SExtractor. The correction has not been included in 
the previous studies, whereas stacking 2-D images with¬ 
out any geometrical corrections could, as is evident, re¬ 
sult in morphologically biased (more centrally concen¬ 
trated) images. We show the comparison of the two me¬ 
dian stacked profiles, which are inclination corrected and 
uncorrected, in Fig. 3b. (For the details of the median 
stacking method, see the following subsection.) We see 
larger discrepancy at the outer radius, where uncorrected 
profiles undergo for ^ 30% of the corrected ones at any 
redshifts, while we see similar profiles at the inner ra¬ 
dius (< 2 kpc). This is because the inclination effect is 
smaller at the rounder inner region. The effects would 
worsen S/N and artificially bias the stacked profiles to¬ 
ward more centrally concentrated ones. Accordingly, in 
what follows we make use of the 1-D profiles with the axis 
ratio and position angle corrections to derive the stacked 
average profiles. 

3.3. Staeking Analysis 

In this subsection we stack the ID light profiles in each 
redshift bin to obtain the composite radial profiles. The 
procedure gives sufficient S/Ns for each pixel, even for 
the outer envelope of the high- 2 ; galaxies. Furthermore, 
the stacking analysis provides us with the radial pro¬ 
files for all filter bands, because an observation field that 
lacks some filter bands (for example, COSMOS field lacks 
F435W, F775W, and F850LP) can be compensated by 
the data taken in other fields with those bands. The full 
coverage of the observed bands would warrant reliable 
SED models. 

The stacking of galaxies is conducted for each filter 
band as follows. We use the inclination corrected 1-D 
imaging data, prepared in the previous section. The 
background objects and contaminants are masked, and 
the masked pixels are not used for the median calcu¬ 
lation when stacking. Since the 1-D conversion relies 
on the center position in the original 2-D image, we use 
galaxy positions listed in the 3D-HST catalog (see third 
and force column in Fig. 2) for the stacking. 

In addition to the local sky subtraction in the previous 
subsection, we compare the sky values of the final stacked 
profile by estimating at different radii (50<r< 75,75 < 
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Fig. 2.— Visual summary of the methodology in Section 3: 1 . sample selection (Section 2.2), 2 . PFS matching (Section 3.1), 3 . sky 
subtraction and 1-D conversion (Section 3.2), 4 . median stacking the 1-D profiles in each redshift bin for each filter band (Section 3.3), and 
5 . Radial SED fit for each pixel (Sections 3A and 4.2). We obtain the radially resolved stellar mass and rest-frame colors for MWs and 
MGs. 


r < 100,100 < r < 125 and 125 < r < 150 pixel), 
and then subtracting the median values from the stacked 
profile, although they are negligible (see appendix for the 
detailed discussion). 

Another concern of stacking analysis is the luminosity 
bias due to redshifts and stellar masses. Since we di¬ 
vide the sample into finite redshift and stellar mass bins, 
the luminosities differ depending on redshifts and stellar 
masses, which would give rise to bias when stacking pro¬ 
files. To avoid the effect, we correct the luminosity for 
each galaxy by multiplying the following constant. 


K{{z), (M*)) 


(1 + ^obs)^ 

(l + (^))3 M*,obs’ 


( 7 ) 


where “obs” represents the observed quantity and the 
brackets are the medians of each redshift bin. It should 
be noted that the correction includes the LC-correction. 
Since the SED fit (implicitly) includes the LC-correction 
when applying the redshift shift, the latter part of Eq. 7 
((M^)/M^^obs) corrects two effects: the LC-correction and 
stellar mass difference. The remaining part of the equa¬ 
tion ((1 -h ^obs)^/(l + (^))^) corrects the luminosity dis¬ 
tance. 

Eigure 4 shows examples of radial S/N profiles for the 


stacked galaxies at ^ ^ 0.7 and ^ ^ 2.8. We find that 
the region with S/N (> 3) reaches at least ^ 10 pixel 
(^ 5 kpc at 2 : ^ 2) from the galaxy center, even in the 
highest redshift bin. The number of stacked samples for 
each redshift bin and filter is summarized in Table 1. 

3.4. SED Fitting 

We calculate radially resolved SEDs for the stacked ra¬ 
dial profiles using EAST (Kriek et al. 2009) and the stel¬ 
lar population models of GALAXEV (Bruzual & Char¬ 
iot 2003), assuming solar metallicity and Chabrier (2003) 
IME. The MW dust attenuation of Cardelli et al. (1989) 
is adopted in the range of 0 < Ay < 4.0 mag by a step of 
0.1. The age is set to range from 0.1 to 10 Gyr, or up to 
the age of the Universe at the galaxy redshift. Redshifts 
are set to the median values of the sample bins. 

The star-formation history is assumed to be an expo¬ 
nentially declining model, SER oc exp(—t/r), where t is 
the time since its star-formation starts and r is the e- 
folding timescale of the SER. It is noted that the ex¬ 
ponential model was first proposed for local elliptical 
galaxies, and recent studies have shown that the de¬ 
layed and truncated models would be consistent with 
observations of high- 2 : star-forming galaxies (Maraston 
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Fig. 3a. — Examples of the PSF profiles (stacked stars provided 
by the 3D-HST team) convolved to match to that of the F160W 
image (red shade). Vertical dashed lines represent the seeing radius 
(~ 0.''09) and physical scale (~ 10 kpc at 2 : ~ 2). The width of the 
red shade in the bottom panel represents an error range within a 
±1% difference from the one of F160W. 




Fig. 3b. — Comparison of the two median light (F160W) profiles 
that are stacked with (solid) and without inclination correction 
(dashed). The profile flux is arbitrarily shifted for each redshift. 
The profiles are similar in the central region (< 2 kpc), but show 
deviation at the outer radius, for ~ 30% at ~ 10 kpc, which would 
bias the profile more centrally concentrated in case of uncorrected 
stacking. Vertical dotted lines represent seeing radius. 




Fig. 4. — Examples of the radial signal-to-noise (S/N) profiles for 
the stacked MWs at {z) = 0.7 (top) and 2.8 (bottom). The vertical 
dashed line represents the PSF radius (FWHM/2) for each redshift. 
The horizontal dashed lines represent the criterion, S/N > 3, for 
radial SED analysis. We conduct the SED fit for each radial pixel 
with S/N > 3. 


et al. 2010; Barro et al. 2013). However, in this study 
we adopt the exponential model because with only the 
photometric data we are not able to discriminate these 
star-formation histories. Instead, we adopt a dispersion 
of the best-fit stellar masses derived with the three star- 
formation histories (exponential, delayed, and truncated 
models) as a typical error of the stellar mass of the ex¬ 
ponential model. The error 0.1 dex) is much larger 
than those (e.g., due to the sky subtraction error). See 
the appendix for the discussion of the uncertainties in 
the stellar mass. It is noted that stellar mass is the 
most reliable parameter to be derived in SED fitting, 
as demonstrated by Wuyts et al. (2012), whereas the 
other parameters wander through the complex parame¬ 
ter space (e.g., well-known age-metallicity-dust degener¬ 
acy; see also Papovich et al. 2001; Shapley et al. 2005). 
The visual summary of the methodology described here 
is shown in Fig. 2. In the next sections, we show the 
results of the radial SED fitting for MWs and MGs. 
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4. RESULTS 

4.1. Integrated SEDs of Staeked Samples 

To study the general properties of the progenitor and 
to assess the accuracy of the SED fitting, we first inte¬ 
grate the stacked radial profile within the central 15 kpc 
radius. By using the stellar population models in Sec¬ 
tion 3.4, we obtain the SEDs of both samples for each 
redshift bin and show them in Eig. 5. We find a good 
agreement for the derived parameters with the median 
values of each sample, shown in Table 2, within error. 
(Eor the calculation of the error, see the following sec¬ 
tion.) 

4.2. Radially Resolved SEDs 

We obtain the radial SED for the stacked samples. Us¬ 
ing the SED templates prepared in Section 3.4, the ra¬ 
dial SED fitting is conducted only for the pixels where 
the fluxes of E850LP, E125W, and E160W have S/N > 3. 
This is because these bands are sensitive to the Balmer 
break at high redshift {z > 1.5). Instead, we verified this 
after the Monte Carlo simulation (see Appendix) and 
found that if the pixel satisfies the criterion, the derived 
parameters have negligible error compared to, for exam¬ 
ple, one caused by the assumption in the SED templates 
(Sec. 3.4). In Eig. 6a, we show the derived surface stellar 
mass density profiles. The stellar surface density is ob¬ 
tained down to ^ 10^ M0kpc“^ at the outer part. The 
error for the stellar mass is estimated from the differ¬ 
ences between the star-formation histories, as described 
in Section 3.4. The error from the photometric uncer¬ 
tainties is comparatively small (^ 10 ^ ATq), and thus we 
neglect them hereafter (see Appendix). 

To conquer the effect of a PSE convolution, we apply 
the method of Szomoru et al. (2010), rather than directly 
deconvolving the images, to the 1-D stellar mass pro¬ 
files. This is because the direct deconvolution could in¬ 
troduce unexpected noise, especially for small size imag¬ 
ing data. Eirst, we derive the best-fit 1-D Sersic pro¬ 
file with GALEIT (Peng et al. 2010), taking account for 
the convolution effect. Then, we add the residual be¬ 
tween the observed and the best-fit profiles to the best-fit 
unconvolved one. Although the residual still remains un¬ 
deconvolved (^ 1/10 of the best-fit model), we verify the 
result by convolving the final “deconvolved” image with 
the original PSE kernel for E160W, and then we find that 
the method reasonably reproduces the observed image, 
except for outermost parts, which hardly affect our final 
results. 

Eigure 6a shows the deconvolved radial profiles for each 
redshift bin, where we see distinct differences between the 
two samples. The stellar mass of MWs is accumulated at 
all radii in a similar way (we mention it as a “self-similar 
way”), from 2; ^ 2.8 to 1.0 with slight surplus increases 
at the outer part (> 2.5 kpc) at 2; < 1. The stellar 
mass of the inner part (r < 2.5 kpc, hereafter the bulge), 
grows by a factor of ^ 7 over the considered redshift 
range. ^ We also show cumulative stellar mass profiles 
in Eig. 6b with non-parametric half-mass-radii. The only 
mild change of half-mass-radius (rh ~ 2.8-5.0 kpc) over 

^ We refer to the inner region (< 2.5 kpc) as “bulge”, because 
the bulge is expected to dominate the region, though for more 
specified definition we need additional analyses, such as the bulge- 
disk decomposition, which is beyond the present scope. 


the redshift supports the self-similar evolution. On the 
other hand, the half-mass radius of MGs evolves from 
^1.9 kpc at 2; > 2 to ^ 4.8 kpc at 2; < 1, by a factor of 
> 2.5. The difference of the stellar mass accumulation 
histories for both samples would become more evident 
when comparing the stellar masses of the bulge and outer 
part in Eig. 6c. The comparable growth of both regions 
for MWs, in the top panel of the figure, supports the self¬ 
similar evolution, while we see in the bottom panel the 
“inside-out” growth for MGs, where the inner massive 
bulge has already appeared at 2: ^ 2.0 (M^^2.5kpc ^ 5 x 
IO^^Mq), after which the mass increases preferentially 
at the outer region. 

Although the evolution trend of is consistent with 
those found in the previous studies, the sizes are slightly 
smaller (for example, re ^ 6 kpc at 2; ^ 0.5 for MGs in 
Patel et al. 2013a, where Ve is the effective radius calcu¬ 
lated by GALEIT). The disagreement would come from 
two possible facts; one from the difference between the 
stellar mass and light profiles, and the other from the dif¬ 
ference of the definition of radius (e.g., non-parametric 
and parametric). The former offset could be originated 
by the effect of color gradients (i.e., higher M/L at the 
central part; e.g., Cappellari et al. 2006), which would 
make a radial profile of stellar mass more centrally con¬ 
centrated. The radial profiles of stellar mass to the H- 
band light ratio for both populations are shown in Eig. 7. 
We find that the slope of M/L varies over the redshift 
range, which suggests that the mass profile is not cor¬ 
rectly reconstructed with only the light profile. We note 
that Szomoru et al. (2013) also found that half-mass-radii 
are ^ 25% smaller than half-light radii for the galax¬ 
ies at 0.5 < 2; < 2.5. Eor the latter, the parametric 
measurement with GALEIT assumes that the outer part 
of galaxies extends to infinity, which would give larger 
half-light radius, while the non-parametric method takes 
account only of really observed stellar mass. The non- 
parametric analysis usually suffers from the fact that the 
outskirts of galaxies at higher redshift could be buried 
under background noise and missed, while the paramet¬ 
ric analysis suffers from the degeneracy of the parameters 
(Sersic index and Ve). We refer the reader to Morishita 
et al. (2014) for the detailed discussion of these issues. 

4.3. Rest-Erame Color Gradients 

In addition to stellar mass profiles, we investigate rest- 
frame UVJ colors. The UVJ colors derived from the 
best-fit SEDs are used to diagnose high-^ quenched and 
star-forming galaxies (Williams et al. 2009; Whitaker et 
al. 2012). We here show the UV J diagrams and U — V 
color profiles for MWs and MGs. 

Eirst, the UV J color diagrams of the stacked ra¬ 
dial profiles are given in Eig. 8. The pixels that en¬ 
ter the top left region, which are bordered with dot¬ 
ted lines, are diagnosed as quiescent, and the others 
are star-forming. The criteria were originally defined 
in Williams et al. (2009) for galaxies at 2; < 2. Mor¬ 
ishita et al. (2014) refined them for the higher redshift 
{z ^ 3) with observed SERs, and we here adopt refined 
ones. (Note that the median specific SERs in Morishita 
et al. (2014) for QGs at 0.5 < ^ < 1.0,1.0 < ^ < 2.0, 
and 2.0 < 2; < 3.0 are 0.03, 0.08, and 0.16 per unit 
Gyr, respectively, while those for star-forming galaxies 
are 0.57, 1.88, and 2.63 per unit Gyr.) The symbols are 
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Fig. 5. — Observed central (< 15 kpc) fluxes (red circles) for MWs (top) and MGs (bottom). The best-fit SED templates calculated 
with FAST are shown with solid lines. 


distinguished according to the distance from the galactic 
center to see the inner (< 2.5 kpc, or bulge) and outer 
(2.5 < r < 10 kpc) colors separately, which is inspired 
by recent findings of the relation between the quiescence 
and compactness of high-z galaxies (Franx et al. 2008; 
Cheung et al. 2012; Woo et al. 2015). We see that at all 
redshift bins the inner pixels stand at the reddest fron¬ 
tier, which is in agreement with the general picture of 
galaxy evolution (Vila-Costas & Edmunds 1992; Wuyts 
et al. 2012; Szomoru et al. 2013). It is noted that al¬ 
though the dust attenuation makes the galaxy color red. 


this does not significantly affect the quiescent sample 
here because dust attenuated galaxies are shifted to just 
right of the quiescent criteria in Fig. 8 (see arrows in 
Fig. 10a; see also Whitaker et al. 2012). 

Figure 9 shows the evolution of radial U — V color, 
which traces the strength of 4000 A and Balmer breaks 
for both populations, where we can conjecture when and 
how the quenching starts. For MWs, we see the bulge 
start to become redder at z ^ 2, compared to the outer 
region. After that, the bulge seems to quench at z ^ 1.0, 
while there is still ongoing star-formation in the outer 
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Major axis (kpc) 



Fig. 6a. — Radial profiles of the stellar surface density of MWs 
(top) and MGs (bottom). The position of each symbol represents 
the half-mass radius, whose value is shown in the caption. The 
vertical dashed lines represent the radius that corresponds to the 
maximum PSF FWHM/2 for the full sample, although the PSF 
convolution effects are corrected (see the text). The bar repre¬ 
sents the typical error for each radius, which is calculated based on 
the difference among the three different star-formation histories. 
Stellar mass within 2.5 kpc from the center (M* 2.5 kpc) total 
stellar mass (M*,tot) for each redshift bin are shown in the caption. 

part, that is, we see the “coexistence” of quenched bulge 
and star-forming disk in MWs. For MGs, on the other 
hand, quenching seems to start at an earlier epoch (z ^ 
2.0), and star-formation activities seem to stop at ^ ^ 
1.0 in whole. We discuss these findings in the following 
section along with the result of the radial stellar mass 
profiles. 

5. DISCUSSION 

In the last two sections we stacked the multi-band 1-D 
light profiles of two populations (MWs and MGs) in each 
redshift bin to obtain the composite radial profiles, and 
then derived the radial SED to see the transitions of the 
radial mass profiles over the redshift range. With the ra¬ 
dially resolved SED profiles, we are now able to discuss 
not only where the stellar mass growth happens, but also 
what mechanisms facilitate the evolution. Eurthermore, 




Fig. 6b. — Cumulative stellar mass profiles of MWs (top) and 
MGs (bottom) for each redshift bin. Each profile is normalized to 
the total stellar mass. The vertical lines and symbols are the same 
as Fig. 6a. The evolution of the half-mass radius is in the inset. 
While the mass profiles of MWs show the self-similar evolution, 
those of MGs show inside-out evolution, where the accumulation 
of stellar mass continues at the outer part of galaxies. 


the stacked profile allows us, for the first time, to evalu¬ 
ate the dispersion of radial profiles of individual galaxies 
around the median profile at a given redshift, as an in¬ 
dicator of morphological variety. Examining this variety 
gives us a clue to the evolution of galaxy morphology, in 
the context of the cosmic star-formation history over the 
redshift and manifestation of the Hubble sequence. 

5.1. Inside-out Mass Growth Versus Self-similar Mass 

Growth 

As shown in Section 4.2, the stellar mass profiles of 
MGs and MWs evolve in different ways. The former 
evolves in an inside-out way, while the latter evolves 
in a self-similar way. In the inside-out scenario, it is 
claimed that massive naked bulges appeared at 2: ^ 2 (or 
even earlier; Nelson et al. 2014), and then it gains the 
rest of mass at the outer part by accreting less massive 
satellite galaxies, as shown in early observation studies 
(e.g., Trujillo et al. 2012). We clearly see in the bottom 
panels of Eigs. 6a and 6b that the massive bulge of the 
MGs (logM^^2.5kpc > 10.5) was formed at 2; ^ 2.4, and 
then the mass growth mostly occurred at the outer part 
(> 2.5 kpc). This can be also found when we compare the 
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Fig. 6c. — Evolution of the stellar mass for MWs (top) and MGs 
(bottom); total mass (black squared), mass within r < 2.5 kpc 
(green circles), and mass in r > 2.5 kpc (yellow triangles). The 
bar represents the typical error of stellar mass at each region. We 
can conjecture two distinct profile evolutions for MWs (self-similar 
way) and MGs (inside-out way). 


mass of the bulge to that of the outer disk in the bottom 
panel of Fig. 6c, where we find that the former does not 
evolve after 2: ^ 2, while the latter continuously evolves 
over the entire redshift range. Furthermore, the observed 
growth of the bulge and total stellar mass suggests that 
more than ^ 75% of the stellar mass is accumulated at 
the outer envelope of galaxies after z ^ 2. Our findings 
for MGs support previous studies with light profiles (van 
Dokkum et al. 2010; Patel et al. 2013a). 

The inside-out, or “two-phase” (Naab et al. 2009; Oser 
et al. 2010), evolution scenario explains the previously 
observed size growth of Tq by a factor of ^ 2-5 (Bezan- 
son et al. 2009; van Dokkum et al. 2010), although the 
scenario is still in dispute (e.g., Carollo et al. 2013; Pog- 
gianti et al. 2013; Riechers et al. 2013). According to the 
scenario, both effective radius and Sersic index are ex¬ 
pected to evolve (Morishita et al. 2014), whereas in this 
work we only measure the non-parametric half-mass ra¬ 
dius to avoid the parametric degeneracy. The formation 
of bulges at z > 2, on the other hand, would be difficult 
to expect in such a simplified and continuous scenario, 
and at the higher redshift violent dynamical mechanisms, 
such as major merger and clump migration, might trig¬ 
ger the formation of bulges (Barro et al. 2014; Dekel et 
al. 2014). 

On the other hand, we find the mass profiles of MWs 
evolving in self-similar way, where the mass accumula¬ 
tion at the inner and outer parts are comparable (see 
the top panel of Fig. 6c). This similarity is a conse¬ 
quence of the continuous bulge growth even at 2; < 1.0, 
while MGs show little or no bulge evolution after 2; ^ 2. 



10 15 20 25 

Major axis (kpc) 



0 5 10 15 20 25 30 35 


Major axis (kpc) 


Fig. 7. — Radial profiles of stellar mass to if-band (F160W) 
light ratio (M*/Lfi6ow) of MWs (top) and MGs (bottom) for 
each redshift. The bar represents the typical error of stellar mass 
at each radius. The vertical lines are same as for Fig. 6a. We see 
variations of M*/Lfi6ow along galaxy radius, which suggest the 
necessity of the radially resolved SEDs to derive the radial stellar 
mass prohles. 


The result is in good agreement with vD13, while there 
are other inconsistent observational evidences for two- or 
three-phase formation scenarios for the MW where the 
bulge and thick disk evolve first, and then the thin disk 
enlarges the galaxy size (Toyouchi & Chiba 2014; see 
also Chiappini 2001). Although we could not elucidate 
the true evolution process only by observing stellar mass 
profiles (because they do not tell us specific dynamical 
processes, e.g., stars radially mixing; Sellwood & Bin- 
ney 2002), merging of galaxies and AGN feedback (Fan 
et al. 2008), in the following give a possible explanation 
for the self-similar evolution, on receiving the results of 
rest-frame colors. 

5.2. Rest-Frame Color Profile as a Diagnosis of 
Quenehing 

In this subsection we investigate the shutdown of star- 
formation activities (quenching) for both populations by 
using the stellar mass and the rest-frame color profiles 
derived in Section 4. The U — V colors (Figs. 8 and 9) 
tell us how the star-forming region becomes quiescent, 
namely, “how the quenching proceeds.” We see that both 
populations quench the star-formation activity from in¬ 
ner to outer radii. The trend is clear in Fig. 9, where the 
inner regions keep reddest at any redshift bins. 
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Fig. 8. — Rest-frame UVJ diagrams for MWs (top) and MGs (bottom). The symbol represents each pixel, and is distinguished with the 
distance from the galactic center (green circles for r < 2.5 kpc and yellow triangles for r > 2.5 kpc). The dotted lines are the boundary 
for quiescence (hatched region) defined in Morishita et al. (2014). The error bars represent the typical error estimated in Section 5.2. The 
evolution of the colors proceeds from the inner to outer regions, which is a common trend for both populations. 
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Fig. 9. — Rest-frame U — V color radial profiles for MWs (top) and MGs (bottom). The bar represents the typical error. The same 
criterion for quiescence as Fig. 8, {U — V) = 1.35, is shown with horizontal dotted lines. 
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For further investigation, we show the color evolution 
of bulges and outer regions of MWs and MGs in Fig. 10a. 
The UVJ criterion is same as in Fig. 8, and the rest- 
frame colors for both samples are derived after summing 
up the pixels in the bulge and outer region. Shown in 
the top panels of Fig. 10b is the evolution oi U — V 
color, where we set the “quenched” region (hatched) as 
{U — V) > 1.35. The error is estimated in the same 
manner for stellar mass (see Section 3.4). The trend is 
similar to what we see in Fig. 8. Although the definition 
of “quench” (or quiescence) is based on both U — V and 
V — J colors, we verify in Fig. 10a that most of the data 
points that have {U — V) > 1.35 are in the quiescent re¬ 
gion of the UV J colors, and therefore we hereafter adopt 
the U — V color as an indicator of quenching. The data 
points of {U — V) > 1.35 mag located outside the qui¬ 
escent hatched region in Fig. 10a, but within the error, 
are shown with open symbols to distinguish them from 
the quenched sample. In the figure, we see the similarity 
of the bulge and outer colors at the highest redshift, for 
each sample. As redshift decreases, on the other hand, 
we find that the bulge becomes redder than the outer 
part (at ^ ^ 1.6 for MWs and 2 : ^ 2.4 for MGs). We 
conjecture two possibilities for bimodality of the colors; 
one is reddening of the bulge through quenching, and an¬ 
other is increase of star-formation activity at the outer 
part, which leads to bluer color. 

To investigate this, we calculate for each epoch the 
expected color at the next epoch by using the best-fit 
SED template and population synthesis model, and then 
compare it with the actually observed color at the next 
epoch. The SED model is prepared with the best-fit pa¬ 
rameters derived in the same manner as the previous sec¬ 
tion, and its evolution is calculated with GALAXEV by 
following the derived SED parameters (SER, r and age) 
at each redshift to the next redshift. Shown in the top 
panel of Eig. 10b, the thin line is the expected evolutions 
of U — V color for each redshift. Then, we obtain the 
model U — V color, {U — V)modeU fhe ith redshift, 
by using the one expected from the {i — l)th redshift to 
compare with the observed U — V color, (1/ —E)o6s, at 
Since GALAXEV does not include quenching functions 
(e.g., gas stripping or AGN feedback), the model and ob¬ 
served quantities should be matched unless such quench¬ 
ing mechanisms work. The comparisons of the colors for 
both populations are depicted in the bottom panels of 
Eig. 10b, with A{U -V) = {U- V)obs -{U- V)model^ 
Erom the figure, we clearly see that the bulges of MWs 
and MGs become unusually red at z ^ 1.6 and 2 : ^ 2.4 
(over a timescale of ^ 1 Gyr), respectively, while the 
outer parts of both populations show neither red nor 
blue excess of color. Therefore, we conclude that the 
bimodality is caused by the shutdown of star-formation 
activities in the inner regions (hereafter referred to as 
“bulge quenching”), while the outer parts evolved as ex¬ 
pected by the best-fit SED models. 

Eor the bulge quenching, we mention several possible 
processes suggested by recent theoretical and observa¬ 
tion studies. The first is the termination of the cold 
gas supply by virial shock at the halo radius. The halo 
masses, Mh, of MWs and MGs at the quenching epoch 
are estimated to be r\j 10^^Mq at 2 ; ~ 1.6 and 
Mh ^ 3 X 10^^ Mq at z ^ 2.4, respectively, based on the 
result of Behroozi et al. (2013). Although the quenched 


halo masses of both populations are not identical, this is 
consistent with, for example, Dekel & Birnboim (2006), 
where cold gas continues to flow into the central part 
of massive halo at higher redshift (Mh > 10 ^^Mq at 
2: ^ 2). At lower redshift, on the other hand, all of it 
is shock-heated, even in the system of lower halo mass 
(Mh ^ 7x 10^^ Mq). The termination of gas flow into the 
central region is also explained by the fact that gas with 
smaller angular momentum preferentially falls into the 
central region while one with larger momentum, which 
accretes at a later epoch, cannot fall and thus remains at 
the outer radius. The second possibility is the one called 
“morphological quenching” (Martig et al. 2009), where 
the formation of the bulge stabilizes the gas kinematics 
and prevents the system from forming stars. Several re¬ 
cent observations with lEUs have found that the bulges 
of high- 2 ; massive galaxies are indeed stabilized with high 
Toomre Q-parameter (e.g., Genzel et al. 2014). The find¬ 
ings of the correlation between the quiescence and mor¬ 
phological compactness also support this scenario (Eranx 
et al. 2008; Gheung et al. 2012; Woo et al. 2015). The 
last possible process is the radiative energy feedback from 
low mass stars in a compact system, which effectively 
prevents the formation of H 2 molecule and stops star- 
formation (Kajisawa et al. 2015). 

In addition, the stellar mass accumulations for both 
populations (Section 4.2) can be investigated with the 
result of rest-frame colors (Section 4.3). As shown in 
Eig. 6c, the bulge of MWs continuously evolves, even 
after its quenching epoch at ^ ^ 1.0. This suggests 
that there should be a mechanism that transports a huge 
amount of stars (^ 5 x 10 ^Mq) into the bulge from the 
outer disk and/or external systems in 1 Gyr or so. The 
bulge of MGs is rapidly formed to be 4 X lO^^M© at 
2: > 2.0, where the SER of MGs peaks (see the follow¬ 
ing section), and after that we see no evolution. The 
outer part of MGs, on the other hand, still continues to 
form stars (^ 5 x lO^^M©) after 2 ; < 1.5, where the star- 
formation seems to stop (Eigs. 8 and 9). The continuous 
growth of the outer envelope is consistently explained 
with ^ 2 merger events of the gas poor satellite galaxy 
with ^ 3 X lO^^M©, which is in good agreement with re¬ 
cent observation studies (e.g., Eerreras et al. 2014). The 
gas poor satellite merger would enlarge the Sersic index 
of galaxy, as was shown in our previous study of massive 
QGs (Morishita et al. 2014). 

5.3. Morphological Variation and Appearance of the 
Hubble Sequence 

We here study the varieties of individual profiles from 
the median stacked profiles and their evolution (see also 
Kajisawa & Yamada 2001). In the top panel of Eig. 11 
we superpose all radial luminosity profiles of galaxies in 
2 ; ^ 0.7 bin, and compare them with the median radial 
profile ^ . The median and individual profiles here are 
the luminosity in E160W-band, since we did not obtain 
the stellar mass profile for each sample. The individual 
light profiles are reduced in the same manner for the 
radial SED, including a mask for background galaxies 
and redshift/stellar mass correction of Eq. (7). 

® The whole samples are available at 
http: //www.astr.tohoku.ac.jp/~mtakahiro/moril5apj/mori_15_rad.pdf 
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Fig. 10a. — Evolutions of rest-frame UVJ colors for MWs (left) and MGs (right). The symbols represent the integrated colors of the 
inner region (r < 2.5 kpc, green circles) and outer region (2.5 < r < 10 kpc, yellow triangles) at each redshift, from the bottom left (z ~ 2.8) 
to the top right (z ~ 0.7). The bar represents the typical error, which is estimated in Section 5.2, for each redshift. The boundary with 
the dashed lines and hatched region are same as for Fig. 8. The effect of dust attenuation (AAy = 1.0 mag) is shown with an arrow. 






Lookback time (Gyr) Lookback time (Gyr) 

Fig. 10b. — Top: evolution of the rest-frame U — V colors of the inner (r < 2.5 kpc; green circles) and outer (2.5 < r < 10 kpc; yellow 
triangles) regions of MWs (left) and MGs (right). The error bar represents the typical error for U — V dX each redshift. The hatched quench 
region is set with (U — V) > 1.35 mag based on Fig. 10a. The data points which locate out of the quiescent hatched region in Fig. 10a 
within the error, whereas (U — V) > 1.35 mag are shown with open symbols to distinguish them from the quenched sample. The model 
predictions ofU — V color (thin lines) are depicted with the best-fit SED profiles at each redshift and the population synthesis model of 
GALAXEV, to see when each region becomes “unusually” red or quenched. Bottom: Gomparison of the observed and model U — V colors 
at each redshift, where we define A(U — V) = (U — V)obs — {U — V)model- The errors are estimated in the same manner for the top panel. 
We see the excess of (U — V)obs in inner regions (r < 2.5 kpc) at 2 : ~ 1.6 and 2 : ~ 2.8 for MWs and MGs, respectively, but no excess 
are seen in the outer regions for either sample. 


We define the variety of the galaxy profile by evaluating 
the distribution of residuals as, 

_ f I /i,median fi, ohs i /q\ 

norm,a: — ^ | 7 5 wJ 

^ Ji,median 

where we set x < 2.5 kpc and 2.5 < x < 10 kpc in the 
following to see the radial dependence of the variety. Sx 
is the total pixel within the range of x. In Fig. 12 we 
show the distribution of Anorm,a: for each redshift bin. 
Then we evaluate the dispersion, a, around the mean of 
the distribution, which is given in each panel. Since the 
scatter of background noise for each light profile could 
affect the quantity, we calculate the error for a by as¬ 


suming the maximum deviation from the best-fit value. 
The results do not change in other filter bands, as long as 
the rest-frame W or longer wavelength bands are used. 

The evolution of a over the redshift range is shown 
in Fig. 13. In the figure, we see that both populations 
have their own characteristic morphological variety. The 
error bars include both of the photometric error and the 
scatter of background noise. For MWs, a of the bulge 
peaks around 2 : ^ 2, and moderately decreases toward 
the lower redshift, whereas at the outer part a remains 
large 0.6). This difference suggests that the bulge of 
MWs has formed at z > 2, while star-formation at the 
outer part has continued rather randomly over the entire 
redshift. This might lead us to an understanding of the 
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Fig. 11.— Examples of the radial light profile for MWs at 2; ~ 0.7. Top: gray scale represents the fraction of radial light (F160W) 
profiles of the galaxies used for stacking. The median profile is shown as a thick solid line with circles. The flux is normalized at the center 
of the median profiles. Middle and bottom panels: Radial light profiles of randomly selected samples (dotted lines), compared with the 
median profiles (solid line; same as in the top panel). The rgb color pictures composited with F814W (blue), F125W (green), and F160W 
(red) bands are shown in each panel with the physical size. The black blot in the bottom rgb image is the masked region. The pictures of 
all samples are available at http://www.astr.tohoku.ac.jp/~mtakahiro/moril5apj/mori_15_rad.pdf. 


self-similar evolution, where the bulge mass grows by a 
factor of ^ 7 from z > 2 to z ^ 0.7 (Section 4.2). To 
interpret the results, some secular process in the galaxies, 
such as stellar migration, would be needed. MGs have, on 
the other hand, the peak at z ^ 2.8 (or higher redshift), 
and the morphological variety rapidly calms down to ^ 
0.2 at 2: ^ 0.7, which is much smaller than those for 
MWs. The finding suggests a very quick transition into 
morphological similarity, in entirety, after the formation 
of the bulge at z ^ 2. 

The results above are well explained in the physical 
context by comparing the evolution of average SFRs with 
those of a in Fig. 13. The median SFRs are estimated 
from the SED fitting for each sample in the 3D-HST cat¬ 
alog, and summarized in Table 2. The SFRs for MWs 
and MGs peak at z ^ 2.0 and 2; ^ 2.8, respectively, 
which are both consistent with the result of the variety 


peaks. It would be reasonably explained by the fact that 
the galaxy morphology begins to have a variety when 
the star-formation activity becomes high at, for example, 
clumpy regions. Although our study is based on the stel¬ 
lar mass, it would not be harsh to speculate that the peak 
redshift also depends on the halo mass of the host galaxy; 
the SFR of the galaxy in massive halos peak at higher 
redshift and rapidly decline, whereas in smaller ones it 
peaks at lower redshift and slowly declines (Behroozi et 
al. 2013; Moster et al. 2013; McDermid et al. 2015). The 
evolution of the morphological variations is consistent 
with the results of the quenching diagnosed by the UV J 
colors (Sections 4.3 and 5.2). The star-formation in disk- 
like (younger) galaxies randomly occurs in both inner 
and outer regions, mainly driven by gas accretion. The 
galaxy merger also increases the morphological variety. 
Although distinguishing which mechanism is the main 
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Fig. 12.— Histograms of the radial variations, Anorm (see the main text for the definition), for MWs (top) and MGs (bottom) of 
r < 2.5 kpc (filled green histograms) and 2.5 < r < 10 kpc (unfilled yellow) for each redshift panel. The dispersion (cr) of distribution, 
which represents the variety of galaxy morphology, is shown in each panel. 


driver of the morphological variety at each epoch would 
be beyond the present study, we note that, according to 
visual inspection of each galaxy image, both mechanisms 
do make the galaxy morphology amorphous. Further ef¬ 
fort is needed to investigate the variety, in conjunction 
with stellar mass and its environment, for example. 

After quenching, there is no such significant star- 
formation because there is no longer cold gas accre¬ 
tion any more (see discussion above), and some physi¬ 


cal mechanisms, such as the dynamical friction of stellar 
components, relax the system in a similar morphology. 
The quenching of galaxies happens at the correspond¬ 
ing redshift of the appearance of morphologically well- 
featured galaxies at 2 ; ^ 1 (Dickinson 2000; Labbe et 
al. 2003; Conselice et al. 2005; Ravindranath et al. 2006). 

Lastly, we stress that the present study does not focus 
on the variety due to any contaminants, such as galaxies 
that do not evolve into MWs at 2 ; ^ 0. Instead, the vari- 
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ety here represents the diversity from “typical” galaxies 
at each redshift epoch selected with the constant num¬ 
ber density, and the time evolution of the morphological 
diversity of the population. 

6 . CONCLUSIONS 

In this study we made use of the HST/WFC3 and 
Advanced Camera for Surveys multi-band imaging data 
taken in CANDELS and 3D-HST survey with total sur¬ 
vey area ^ 900 arcmin^ and 5cr source magnitude H ^ 
27.0 ABmag, which allowed us to sample galaxies with 
> 2 X 10^ Mq at 0.5 < 2 ; < 3.0. Based on the con¬ 
stant cumulative number density method, we extracted 
the progenitors of the Milky Way (MWs) and massive 
galaxies (MGs) at the redshift range by following the stel¬ 
lar mass criteria, log(M*/M 0 ) = 10.66 — 0.045z — 0.13z^ 
and log(M^/MQ) = 11.19 — O.O 682 ; — 0.04z^, respectively 
(Fig. 1). Then, we stacked them in each redshift bin 
to obtain the average light and color profiles with full 
care for the sky subtraction and the contamination of 
the neighboring objects (Fig. 2 ). Furthermore, we inves¬ 
tigated the noise, systematic uncertainties (e.g., offsets of 
sky subtraction), and correction of axis ratio, which had 
not been included in the previous studies. After reduc¬ 
ing such noise and bias, we derived the radially resolved 
SEDs (radial SEDs) for the stacked profiles with suffi¬ 
cient S/Ns to obtain the radial profiles for stellar mass 
and rest-frame colors. Based on the present results, we 
discussed not only the general properties and the mor¬ 
phological diversity for both populations, but also possi¬ 
ble mechanisms for cessation of star-formation activities 
(“quenching”) and stellar mass growth at later epoch. 
Our conclusions are as follows; 

1 . MWs accumulated the stellar mass in a self-similar 
way in the bulge and disk. After its quenching, the bulge 
contributed to grow by ^ 5 x 10 ^ Mq after z ^ 1. 

2 . On the other hand, MGs accumulated the stellar 
mass in an inside-out way, obtaining more than 75% of 
the total stellar mass at the outer part (> 2.5 kpc) after 
the rapid formation phase of the massive bulge (M^ ^ 
4 X IO^^Mq) at ^ > 2 (Figs. 6 a, 6 b and 6 c). 

3. Quenching was found to occur strongly, depend¬ 
ing on the stellar (and halo) masses of galaxies and 
bulges. The finding suggests the evidence of bulge- 
related quenching mechanisms, such as morphological 
quenching and termination of gas flow by shock heating 
(Figs. 10a and 10b). 

4. By comparing the median and individual light pro¬ 
files, we evaluated the evolution of the varieties of galaxy 
morphology for the first time. The varieties are rele¬ 
vant to the observed star-formation activities (SFRs and 
rest-frame colors) and the appearance of morphologically 
well-featured galaxies (Fig. 13). 

In this study, we investigated the evolution of the gen¬ 
eral properties and diversity for the two populations, the 
MW-like and massive galaxies, in view of their progeni¬ 
tors as a function of the age of the universe. The popula¬ 
tions have the different epoch of bulge quenching, inside- 
out and self-similar growth of the stellar mass, similarity 
and variety of the morphology, and manifestation of the 
two main structures (bulge and disk) while undermining 
star-formation. The corroboration for the origins giving 
rise to such dichotomy would provide us clues to further 


understanding of galaxy evolution in the unified frame¬ 
work. 
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APPENDIX A 

SYSTEMATIG ERROR GAUSED BY SKY 
SUBTRAGTION 

Since SED fitting derives the best-fit parameters based 
on spectral colors, the wrong sky subtraction on each 
filter image would give significantly biased results. Al¬ 
though space based imaging data, compared to the 
ground-based one, rarely suffers from the background 
light, the zodiacal light still seems to have a nonnegligible 
effect. To confirm whether the sky background value is 
appropriately estimated, we conduct the following test. 
First, we estimate the sky background value at different 
radii, 75, 100, and 125, with a thickness of 25, in unit 
of pixel, for stacked 1-D profiles. (For details about the 
sky subtraction for each unstacked sample, please see the 
main text.) Then, we compare the estimated those sky 
values, including no-sky subtraction (original stacked im¬ 
age), and adopt the median value as the best sky value 
(skybest) and the farthest value from skybest as the max¬ 
imum error (Asky). 

We show the MWs and MGs light profiles (F160W) in 
Figs. A1 and A2, respectively. In the top panels, we show 
the radial light profiles with error bars. The error bars 
include both of data error (Eq. 4) and offset originated 
from Asky. The bottom panels show the cumulative 
light profiles of the best sky subtracted (solid) and those 
with shifts for Asky (dashed). The differences between 
the radial profiles at < 20 kpc are small, even for the 
samples at 2 : > 2, where the sky background noise would 
be dominant. Since our analysis of the radial SED uses 
the inner parts (< 20 kpc at 2 : < 2 and < 10 kpc at 
2 ; > 2) because of the S/N limit, we conclude that our 
sky estimation is robust and the variation of the values 
does not change the final results at all. 

APPENDIX B 

ERROR ESTIMATION IN STELLAR MASS 

We here estimate the error in the stellar mass derived 
by radial SED with the criterion of S/N> 3 for F814W, 
F125W and F160W. Since the SED fit includes several 
parameters that would be degenerate, we estimate the 
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Fig. 13.— Evolutions of morphological variety, cr, of the r < 2.5 kpc (green solid line with circles) and 2.5 < r < 10 kpc (yellow dashed 
line with triangles) regions along with redshift for MWs (top) and MGs (bottom). The error for cr, which is calculated based on the 
photometric error in F160W image, is negligible and smaller than the symbol size. The sample in the shade region (2.5 < 2 : < 3.0) of the 
top panel (MWs) has weaker completeness (but still > 75%) and might be biased toward more luminous galaxies (and higher cr), although 
excluding the data point does not change the trend. Average star-formation rates (gray squares with dotted lines) are derived from the 
3D-HST catalog. The error bar for the SFR represents its dispersion. 
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Fig. Al.— Top: radial light (F160W) profiles of the MWs. The 
vertical axis is the observed flux, which is not corrected for the 
redshift effect. The error bars are calculated by combining Eq. (4) 
in the main text and Asky. Bottom: cumulative light profiles of 
the same sample, for 2 : ~ 0.7-2.8 from top to bottom. Profiles 
are subtracted with the median sky values in the caption and are 
shown as solid lines, while those subtracted and added the largest 
sky values are dashed lines. 




errors with the Monte Carlo test. In the test, we ran¬ 
domly generate the error based on the gaussian distribu¬ 
tion with cr, which is estimated as Eq. (4), for each filter 
band. Then, we shift the observed fluxes for the errors 
above, and derive the best-fit SED for them. We repeat 
this 500 times for each pixel, and estimate the root mean 
square of the best-fit parameter, which we define as the 
error of the parameters. Examples of 10 profiles for MWs 
at {z) = 0.7 and {z) = 2.8 are shown in Eig. Bl, although 
we do not correct the PSE effect for them. Using the cri¬ 
terion of S/Ns > 3, we find that the error is negligibly 
small (< 3 X 10^ Mq at z ^ 0.7 and < 3 x IO^Mq at 
2 ; ^ 2.8), and hardly changes our final results. There¬ 
fore, in the main text we adopt the scatters calculated 
from the different star-formation histories (exponential 
decline, delayed, and truncated) as the typical errors for 
stellar mass and rest-frame colors. 

























































Milky-Way-Like and Massive-Galaxy Progenitors at 0.5 < z < 3.0 


19 




-- log^*.2.5kpc=9.39, 

, logM^_ 

tot =9-59 

-- logAr,_2.5^p,=9.4, 


,t=9.61 

-- log^*.2.5kpc=9.4, 


,t=9.61 

-- logA^*.2.5kpc=9.4, 

logM^^^ 

,t=9.61 

-- log^*.2.5kpc=9.41, 

, logM^_ 

tot =9-6 

-- logA^*.2.5kpc=9.4, 

logM^^^ 

,t=9-6 

-- log^*.2.5kpc=9.4, 


,t=9.61 

-- logA^*.2.5kpc=9.4, 

logM^^^ 

,t=9.61 

-- l0g^*,2.5kpc=9.4, 


,t=9-6 

-- log^*.2.5kpc=9-4, 

logM^^^ 

,t=9.6 

- Median 




10 15 20 

Major axis (kpc) 


Fig. B1.— Examples of 10 sets of stellar mass profiles (dotted 
lines) derived after Monte Carlo tests for MWs at 2 ; ~ 0.7 (top) and 
2 ; 2.8 (bottom). The solid lines with error bars are the median 

profile and the median absolute deviation at each radius. The PSF 
effect is not corrected. 
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TABLE 3 

SExtractor parameters for the object 

MASK AND SKY ESTIMATION. 


Parameter 

Value 

Unit 

DETEGT_MINAREA 

20 

pixel 

DETEGT_THRESH 

1 

(7 

DEBLEND_NTHRESH 

1 


DEBLEND_MINGONT 

1 


BAGK_SIZE 

250 

pixel 

BAGK_FILTERSIZE 

1,1 


BAGKPHOTO_TYPE 

GLOBAL 


BAGKPHOTO.THIGK 

5 

pixel 
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